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ABSTRACT 

In this paper we report the first observation in the Southern hemisphere of an energy dependence 
in the Galactic cosmic ray anisotropy up to a few hundred TeV. This measurement was performed 
using cosmic ray induced muons recorded by the partially deployed IceCube observatory between 
May 2009 and May 2010. The data include a total of 33xl0 9 muon events with a median angular 
resolution of ~ 3° degrees. A sky map of the relative intensity in arrival direction over the Southern 
celestial sky is presented for cosmic ray median energies of 20 and 400 TeV. The same large-scale 
anisotropy observed at median energies around 20 TeV is not present at 400 TeV. Instead, the high 
energy skymap shows a different anisotropy structure including a deficit with a post-trial significance 
of -6.3 <t. This anisotropy reveals a new feature of the Galactic cosmic ray distribution, which must 
be incorporated into theories of the origin and propagation of cosmic rays. 

Subject headings : cosmic rays — anisotropy 
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1. INTRODUCTION 

During the last decades, Galactic cosmic rays have 
been found to have a small but measurable energy depen- 
dent sidereal anisotropy in their arrival direction distri- 
bution with a relative amplitude of order of 10 -4 to 10 “ 3 . 
The first comprehensive observation of the cosmic ray 
sidereal anisotropy was provided by a network of muon 
detectors sensitive to cosmic rays between 10 and several 
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hundred GeV (Nagashima et al. 1998). More recent un- 
derground and surface array experiments in the North- 
ern hemisphere have shown that a sidereal anisotropy 
is present in the TeV energy range (Tibet Air Shower 
gamma (AS 7 ) array (Amenomori et al. 2006), Super- 
Kamiokande (Guillian et al. 2007), Milagro (Abdo et al. 
2009), and ARGO-YBJ (Zhang 2009)). Furthermore, the 
IceCube neutrino observatory reported the first observa- 
tion of a cosmic ray anisotropy in the Southern sky at 
energies in excess of about 10 TeV (Abbasi et al. 2010a). 
The cosmic ray anisotropies reported by IceCube showed 
that the large scale features appeared to be a continua- 
tion of those observed in the Northern hemisphere. 

At high energies, the Tibet AS 7 collaboration re- 
ported an observation for primary energies ~300 TeV 
to be consistent with cosmic ray isotropic inten- 
sity (Amenomori et al. 2006), while the EAS-TOP col- 
laboration reported a sharp increase in the anisotropy for 
primary energies ~370 TeV (Aglietta et al. 2009). At the 
time of the writing of this paper the observations in the 
Northern hemisphere do not provide a coherent global 
picture of the sidereal anisotropy at high energy. 

The origin of the anisotropic distribution in the arrival 
direction of Galactic cosmic rays over the entire celestial 
sky is still unknown. If there is a relative motion of the 
solar system with respect to the cosmic ray plasma, then 
this would produce a well defined anisotropy. For exam- 
ple, if cosmic rays are at rest with respect to the galactic 
center, a dipole anisotropy would be expected. The mag- 
nitude of the anisotropy is calculated to be 0.35% with an 
apparent excess of cosmic ray counts toward the direction 
of solar Galactic rotation (a = 315° ,<5 =48°) and a deficit 
in the opposite direction (a = 135°, <5 =-48°). Such a 
dipole anisotropy is referred to as the Compton-Getting 
effect (Compton &; Getting 1935). Neither the ampli- 
tude nor the phase expected from the Compton-Getting 
effect are consistent with the cosmic ray anisotropy 
observations (IceCube (Abbasi et al. 2010a), Tibet Air 
Shower gamma (AS 7 ) array (Amenomori et al. 2006), 
Milagro (Abdo et al. 2009)). Moreover, the observed 
sidereal anisotropy is not consistent with a simple 
dipole (Abbasi et al. 2010a). It is worth noting that since 
the reference frame of the Galactic cosmic rays is not 
known, it is reasonable to assume that the Compton- 
Getting effect could be (at most) one of several contri- 
butions to the cosmic ray anisotropy. 

While the origin of the anisotropy is not understood, 
it has been speculated that it might be a natural conse- 
quence of the distribution of cosmic ray Galactic sources, 
in particular nearby and recent supernova remnants 
(SNR). The discreteness of such sources, along with cos- 
mic ray propagation through a highly heterogeneous in- 
terstellar medium, might lead to significant fluctuations 
of their intensity in space and time and, therefore, . to 
an anisotropy in the arrival direction of cosmic rays at 
Earth (Erlykin & Wolfendale 2006). This speculation is 
challenged by Butt (2009), who points out that the ob- 
served anisotropy is of low intensity, whereas the high 
energy cosmic rays from such sources would escape the 
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galaxy relatively quickly, leading to high anisotropy. 

The study of the cosmic ray arrival distribution 
might provide hints into the properties of cosmic 
ray propagation in the turbulent interstellar magnetic 
field (Beresnyak et al. 2011). While at TeV energies 
it is speculated that propagation effects could either 
generate large scale anomalies in their arrival direc- 
tion (Battaner et al. 2009) or produce localized excess 
regions (Maikov et al. 2010), depending on the turbu- 
lence scale and diffusion properties, it is still not clear 
whether such models would be able to explain the obser- 
vations at higher energies. 

In this paper we present the analysis of cosmic 
ray data collected by the IceCube observatory, which 
we use to extend the observations of the Galactic 
cosmic ray anisotropies by IceCube (Abbasi et al. 
2010a), (Abbasi et al. 2011) up to several hundred TeV. 
The analysis procedure is described in Section 2 and the 
anisotropy in sidereal reference frame is shown in Sec- 
tion 3. Section 3 describes an experimental procedure 
to verify that the observed sidereal anisotropy is not an 
artifact of the analysis procedure, using the arrival dis- 
tribution of cosmic rays as a function of the angular dis- 
tance from the Sun. In this coordinate system, a dipole 
effect is expected such that the cosmic ray count rate 
is higher toward the direction of Earth’s motion around 
the Sun and lower in the opposite direction. The ex- 
perimental systematic uncertainties on the anisotropy in 
sidereal coordinates are described in Section 4 and the 
conclusions are summarized in Section 5. 

2. ANALYSIS 

2.1. Data and Reconstruction 

IceCube is a neutrino observatorj- located at the geo- 
graphic South Pole. During the 2009-2010 austral sum- 
mer, the partially deployed detector was equipped with 
3,540 Digital Optical Modules (DOMs) buried between 
about 1.5 and 2.5 km below the surface of the ice along 
59 vertical strings (Abbasi et al. 2009). The IceCube 
physics runs in the 59-string configuration (IceCube-59) 
started on May 20, 2009, and ended on May 30, 2010. 
IceCube observes relativistic charged particles by detect- 
ing the Cherenkov light produced as they travel through 
the ice. In particular the observatory is sensitive to the 
charged particles produced by neutrino interactions in- 
side the ice, as well as the muons created in the cosmic 
ray air showers. 

In order to reject the random signals derived from the 
~500 Hz dark noise rate from each DOM, a local coin- 
cidence was required between neighboring DOMs with a 
coincidence time interval of ± 1,000 ns. A trigger was 
then produced when eight or more DOMs in local coinci- 
dence detected photons within 5,000 ns. The trigger rate 
in IceCube-59, predominantly from muons produced in 
cosmic ray air showers, ranged from a minimum of about 
1,600 Hz in the austral winter to a maximum of about 
1,900 Hz in the austral summer. This modulation is due 
to the large seasonal variation of the stratospheric tem- 
perature, and consequently the density, which affects the 
decay rate of mesons into muons (Tilav et al. 2010). 

All recorded events were processed using a coarse on- 
line fit to their trajectories (Ahrens et al. 2004). To re- 
fine the directional estimate, the coarse fit was used to 


seed an online likelihood-based reconstruction, which was 
applied if ten or more optical sensors were triggered by 
the event. The average rate of the events that passed the 
likelihood-based reconstruction ranged from a minimum 
of about 1,150 Hz to a maximum of about 1,350 Hz. All 
the events collected and processed by the IceCube obser- 
vatory were stored in a compact Data Storage and Trans- 
fer format, or DST, and shipped North through satellite 
link (see Abbasi et al. (2011) for details). This analysis 
uses all events with likelihood directional reconstruction 
stored in the DST data format, collected within one full 
calendar year from the beginning of the run on May 20, 
2009. After rejecting short data runs we ended up with 
33 x 10 9 events, corresponding to a detector livetime of 
324.8 days. The events have a median angular resolution 
of about 3° and a median energy of the cosmic ray par- 
ent particles of about 20 TeV. It is worth noting that this 
angular resolution is a property of this data sample and 
the applied reconstruction algorithms; reduced data sam- 
ples using more advanced reconstructions for high energy 
neutrino searches have a typical angular error less than 
1° (Abbasi et al. 2010b). 

To measure an anisotropy of order 10 -4 to 10~ 2 3 , it is 
necessary to eliminate any background effects that could 
mimic such an observation. Due to its unique location 
at the geographic South Pole, the IceCube observatory 
has full coverage of the Southern sky at any time of the 
year. Therefore, seasonal variations in the muon inten- 
sity occur uniformly across the entire field of view and 
do not affect the local arrival direction distribution of the 
reconstructed events (Abbasi et al. 2010a). The main ef- 
fect that needs to be accounted for is due to the geomet- 
rical shape of IceCube: the hexagonal geometric struc- 
ture of the observatory introduces a strong asymmetry 
in the local azimuth distribution of events (Figure 1). 
Non-uniform time coverage caused by detector downtime 
and run selection reduces the total detector livetime by 
about ~ 10%, preventing the complete averaging of the 
local coordinate asymmetry over one year and generat- 
ing spurious variations in the arrival directions of cosmic 
rays in celestial coordinates. To remove this effect, the 
asymmetry in the local azimuthal acceptance (shown in 
Figure 1-b) is corrected by re-weighting the number of 
events from a local azimuth bin to the average number 
of events over the full range of the local azimuth distri- 
bution. This re- weighting is applied in four zenith bands 
with approximately the same number of events per band 
due to the detector azimuth distribution variation with 
zenith angle (Abbasi et al. 2010a). 

2.2. Estimation of Cosmic Ray Energy 

Since IceCube detects cosmic rays indirectly through 
the observation of muons produced in extensive air show- 
ers, the energy of the cosmic ray primary particle is 
estimated based on the total amount of light seen by 
the detector, which is a function of the number and en- 
ergy of detected muons. Muons produced in the atmo- 
sphere propagate through the ice losing energy via ion- 
ization and stochastic processes such as pair production, 
bremsstrahlung and photo- nuclear interactions. The sec- 
ondary charged particles produced by these processes 
emit Cherenkov light. The number of emitted photons is 
proportional to the total energy of the secondaries. By- 
detecting photons, it is possible to estimate the energy 
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Fig. 2. — The average logarithm of the cosmic ray primary energy 
as a function of N& and Zenith angle, as obtained from simulation. 
The Y-axis is the logio of AT ch , the X-axis is the cosine of the 
reconstructed Zenith angle of the event while the color scale is the 
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Fig. 1. — Figure (a) shows the complete IceCube 86-string con- 
figuration, Circles filled in blue represent the IceCube 59-string 
configuration which is the main configuration used in this paper. 
Figure (b) shows the azimuth distribution for the whole data set. 
It shows the number of events vs. the azimuth of the arrival direc- 
tion of the primary cosmic ray particle. The horizontal red line is 
the average number of events for the distribution. 

lost by the muons and therefore the muon’s energy within 
the volume instrumented with optical sensors. However, 
the total energy of the detected muons is only a fraction 
of the original cosmic ray primary energy, while the rest 
is mostly dispersed into the electromagnetic component 
of the air shower. As a consequence, the natural fluc- 
tuations that arise in the development of the extensive 
air showers limit the resolution of the estimate of the 
primary energy that one can make using muons in ice. 

The uncertainty in the cosmic ray energy estima- 
tion has been modeled with a full simulation of cos- 
mic rav interactions in the atmosphere using COR- 
SIKA (CORSIKA 2009) with SIBYLL hadronic interac- 
tion model (Version 2.1) (Engel 1999) together with the 
composition and the spectrum of primary cosmic rays as 
described in Horandel (2003). Muons were propagated 
through the ice with the Muon Monte Carlo (MMC) 


mean of the logarithm of the cosmic ray primary energy for each 
bin obtained from simulation in GeV. The first energy band with 
median energy of 20 TeV is all the events selected belov the dashed 
line, while the second energy band with median energy of 400 TeV 
contains events selected between the continuous lines. 

propagator (Chirkin & Rhode 2004), and a full detector 
simulation was performed on those events. 



FIG. 3. — The fraction of events vs. the logarithm of primary 
energy (in GeV) for the two selected energy samples (see text). 
The low energy sample contains events with a median energy of 
20 TeV (squares) and the high energy sample contains events with 
a median energy of 400 TeV (triangles). The energy distributions 
were determined using a full simulation of cosmic ray interactions 
in the atmosphere, of muons propagation through the ice and of 
the IceCube-59 detector. 

In this analysis the estimate of the cosmic ray energy 
is based on the number of DOMs hit by Cherenkov pho- 
tons (i.e. number of channels, or N c h). The downward 
muons reaching IceCube with a large zenith angle 0 have 
to cross a larger slant depth than vertically propagat- 
ing muons, and so the set of horizontal events naturally 
excludes lower primary energy cosmic rays. This intro- 
duces a zenith angle dependence of the relation between 
ATch and the primary particle energy. Therefore, a two- 
dimensional cut in N c h and 6 is used. Figure 2 shows the 
distribution from simulation of the cosmic ray primary 
particle energy with respect to N c h as a function of cos#. 
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Fig. 4. — The number of events seen by IceCube vs. the logarithm 
of primary energy (in GeV) using the composition model described 
in Horandel (2003). Fractional contributions of proton, helium, and 
iron are shown as well. At 20 TeV, the spectrum is dominated by 
the proton fractional contribution of ~ 70%, while at 400 TeV that 
fraction will have decreased to ~ 30%. The energy distributions 
were determined using a full simulation of cosmic ray interactions 
in the atmosphere as described in this section. 

The figure shows that for a given range of iV C h, vertical 
events (i.e. cos 0 « 1) are dominated by cosmic rays with 
lower average energy than horizontal events (i.e. cos0 
« 0.3) due to the larger ice thickness the muons would 
go through before triggering the detector. We identified 
regions of constant primary energy in {N c h, cos0), de- 
limited with the black lines in Figure 2, in order to select 
two event samples at energies with minimal overlap and, 
at the same time, with the maximum possible number of 
events in the high energy sample. The low energy sam- 
ple was obtained by selecting all events below the dashed 
line in Figure 2, and the high energy sample by selecting 
events between the solid lines in the figure. 

Figure 3 shows the simulated primar 3 T energy distri- 
butions for the two event samples. The estimate of the 
primary cosmic ray energy has a resolution of about 0.5 
in the logarithmic scale. The uncertainty of the primary 
energy estimate is dominated by the fluctuations in the 
air showers. The low energy sample over the Southern 
sky contains 21 x 10 9 events; assuming the composition 
described by Horandel (2003) and shown in Figure 4. 
The median primary particle energy of the low energy 
sample is 20 TeV, with 68% of the events are between 
4 — 63 TeV. The high energy sample contains 0.58 x 10 9 
events. The median primary particle energy of the high 
energy sample is 400 TeV, with 68% of the events are 
between 100— 1, 258 TeV. 


Ii = 


Ni(a y 6) 

(Ni(5)) a ' 


( 1 ) 


where U is the relative intensity for each bin of angu- 
lar equatorial coordinates (a, 5), Ni is the number of 
events in bin i, and ( Ni ) is the average number of events 
for the bins along the same iso-latitude as bin i (with 
the same declination 5). The sky maps in this analy- 
sis are produced using the Hierarchical Equal Area Iso- 
Latitude Pixelization (HEALPix) libraries (Gorski et al. 
2005). HEALPix subdivides the unit sphere into quadri- 
lateral pixels of equal area. In this analysis, the maps 
contain pixels that correspond to an angular resolution 
of- 3°, which approximately corresponds to the angular 
resolution of the detector. 

Figure 5 show the maps of the relative intensity in 
cosmic ray arrival direction in sidereal reference frame 
(equatorial coordinates), for the low and high energy 
samples, respectively. The color scale in the figures rep- 
resents the relative intensity as described in eq. 1. The 
observed sidereal anisotropy appears to evolve as a func- 
tion of energy and the anisotropy pattern observed at 400 
TeV shows substantia] differences with respect to that 
observed at 20 TeV. Note that in the maps only the pixels 
below declination angle of -25° are shown. Pixels above 
declination of -25° are masked due to the degradation 
of the angular resolution at higher declinations. Such 
degradation is to be expected because of the poorer sta- 
tistical power and the domination by mis-reconstructed 
events (Abbasi et al. 2011). 



3. RESULTS 

3.1. Sidereal Anisotropy 

In order to investigate the cosmic ray arrival direction 
distribution, we determine the map of deviation from 
isotropy by calculating the relative intensity distribution 
after azimuthal re- weighting of the arrival directions of 
the data as described in the previous section. The cos- 
mic ray arrival direction distribution is dominated by the 
zenith angle dependence of the muon flux. The zenith an- 
gle dependence is a result of a varying overburden for the 
muons through the ice. Therefore, the flux for each bin 
is normalized within each zenith band (or, equivalently 
at the South Pole, each declination band): 


Fig. 5. — Figure (a) shows the IceCube cosmic ray map of the first 
energy band (median energy of 20 TeV) for the relative intensity 
in right ascension a. Figure (b) shows the IceCube cosmic ray 
map for the second energy band (median energy of 400 TeV) of 
the relative intensity in right ascension a. 

In order to characterize quantitatively the general 
structure of the anisotropy, we proceed as follows. For 
each row of pixels in the map, a 24-bin histogram is made 
from the relative intensity values of the pixels (where 
each pixel’s value is included in the bin which contains 
the right ascension of the center of the pixel). The rows 
are spaced approximately every —3 degrees in declina- 
tion, and the histograms are constructed down to decli- 
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nation -72 degrees (beyond which the number of pixels 
per declination band is less than the number of bins in 
the histogram) . The binned relative intensity data were 
then fitted to a harmonic function of the form 
2 

22 A j C03 \J ( a -<t>j)] + B, (2) 

3 = 1 

where j is the harmonic term order (i.e. dipole for j=l, 
quadrupole for j= 2), Aj is the amplitude of the j th har- 
monic term, <j)j is the phase of the j th harmonic term, a 
is the right ascension, and B is a constant. The results of 
this fit are shown in Table 1 and Table 2 for the low and 
high energy samples, respectively. In addition, in order 
to quantify the sidereal anisotropy over the whole South- 
ern hemisphere, the anisotropy profile in right ascension 
is measured by accumulating the relative intensity dis- 
tribution from the declination belts. The error bars were 
obtained by propagating the statistical errors from each 
declination belt. Figure 6 show the projections in right 
ascension of the cosmic ray relative intensity in sidereal 
reference frame, for the low and high energy samples, re- 
spectively. The lines in the figures represent the fit to 
the first and second harmonic terms of eq. 2, and the fit 
results are shown in Table 3 together with the x 2 /ndof 
values for the first and second harmonic fits, in addition 
to the number of events used in the right ascension pro- 
jections. While the x 2 /ndof indicates that the fits do not 
completely describe the data, we found that even fitting 
up to the sixth harmonic does not completely fit all of 
the structures, so we use here only the first and second 
harmonics as a general characterization of the anisotropy. 


TABLE 1 

Harmonic fit values per declination band for the energy 

BAND CENTERED AT 20 TeV. 


Dec. 

Mean 

Ai ± (stat.) 

io- 4 

^>i± (stat.) 

f°l 

± (stat.) 
IO -4 

02 ± (stat.) 

L°] 

-24 

7.1 ± 1.0 

37.3 ± 8.1 

3.2 ± 1.0 

303.5 ± 9.0 

-27 

8.4 ± 0.9 

35.6 db 6.0 

2.1 ± 0.9 

321.3 db 11.8 

-30 

8.7 ± 0.7 

45.4 ± 4.7 

4.0 ± 0.7 

306.6 ± 5.1 

-33 

8.6 ± 0.7 

50.5 ± 4.3 

3.6 ± 0.7 

294.6 ± 5.0 

-36 

9.3 ± 0.5 

51.2 ± 3.3 

3.1 ± 0.5 

299.1 ± 5.0 

-39 

8.3 ± 0.5 

52.9 ± 3.4 

2.1 db 0.5 

299.6 ± 6.6 

-42 

9.6 ± 0.4 

51.1 ± 2.6 

3.1 ± 0.4 

301.8 ± 4.0 

-45 

9.3 ± 0.4 

57.4 ± 2.8 

3.0 ± 0.5 

305.9 ± 4.2 

-48 

8.0 ± 0.4 

56.7 ± 2.8 

2.7 ± 0.4 

304.3 ± 4.0 

-51 

7.9 ± 0.4 

57.2 ± 2.8 

2.5 ± 0.4 

293.0 ± 4.3 

-54 

8.0 ± 0.4 

55.9 ± 2.6 

2.3 ± 0.4 

297.9 ± 4.5 

-57 

7.9 ± 0.4 

60.8 ± 2.7 

1.8 ± 0.4 

303.3 ± 5.6 

-60 

7.9 ± 0.4 

52.7 ± 2.6 

2.0 ± 0.4 

300.4 ± 5.3 

-63 

7.7 ± 0.4 

49.9 ± 3.3 

1.8 ± 0.4 

307.1 ± 6.7 

-66 

7.3 ± 0.4 

51.0 ± 2.9 

4.1 ± 0.4 

293.2 ± 2.7 

-69 

5.7 ± 0.4 

50.8 ± 4.2 

4.9 ± 0.4 

282.4 ± 2.4 

-72 

5.7 ± 0.4 

38.8 ± 4.0 

3.6 ± 0.4 

301.7 ± 3.2 


Note. — First and second harmonic fit values per declination 
for the first energy band. 


3.1.1. Significance 

Figure 7-a shows the significance map for the 20 TeV 
energy, while Figure 7-b shows the significance map for 
the 400 TeV energy. The significance skymaps are cal- 
culated using the direct integration method with a time 
window of 24 hours and an optimized smoothing as de- 
scribed in Abbasi et al. (2011). The smoothing is then 


TABLE 2 

Harmonic fit values per declination band for the energy 
BAND CENTERED AT 400 TeV. 


Dec. 

Ai ± (stat.) 

<f> i± (stat.) 

A 2 d: (stat.) 

<^2d: (stat.) 

Mean 

10" 4 

i°L, 


10- 4 

[°l a 


-24 

9.6 ± 3.1 

248.1 

± 18.6 

5.4 i 3.1 

143.6 

TT6.6 

-27 

1.1 ± 3.0 

245.7 

± 15.8 

6.5 ± 3.0 

158.1 

± 13.2 

-30 

5.1 ± 2.6 

238.9 

± 29.6 

3.0 ± 2.6 

146.9 

± 25.2 

-33 

3.9 ± 2.7 

255.9 

± 37.8 

2.0 ± 2.6 

205.3 

± 37.6 

-36 

9.6 ± 2.4 

217.0 

d= 14.2 

6.2 ± 2.4 

171.5 

± 10.9 

-39 

9.5 ± 2.4 

246.9 

± 14.3 

6.5 ± 2.4 

144.2 

± 10.5 

-39 

9.5 ± 2.4 

246.9 

± 14.3 

6.5 ± 2.4 

234.2 

± 10.5 

-42 

4.2 ± 2.2 

246.2 

± 30.1 

2.5 ± 2.2 

231.3 

± 25.4 

-45 

1.2 ± 2.5 

311.4 

± 115.6 

2.8 ± 2.5 

110.4 

± 25.1 

-48 

1.4 ± 2.3 

181.0 

± 95.6 

3.6 ± 2.3 

154.2 

± 18.2 

-51 

3.7 ± 2.4 

236.7 

d= 38.2 

2.0 ± 2.4 

156.8 

± 35.6 

-54 

5.5 ± 2.4 

220.8 

± 25.8 

1.5 ± 2.5 

142.5 

± 46.8 

-57 

1.4 ± 2.6 

228.8 

db 112.1 

3.7 ± 2.6 

165.0 

± 21.9 

-60 

3.9 ± 2.6 

359.8 

± 38.5 

7.4 ± 2.6 

161.0 

± 10.2 

-63 

2.6 ± 3.4 

13.0 db 72.8 

3.2 ± 3.3 

148.6 

± 29.6 

-66 

1.3 ± 2.9 

143.4 ± 127.8 

5.3 ± 3.0 

107.5 

± 15.9 

-69 

1.0 ± 3.4 

304.5 db 188.2 

4.2 ± 3.4 

227.9 

± 23.2 

-72 

6.8 ± 3.4 

174.8 ± 28.4 

6.7 ± 3.4 

152.5 

± 14.5 


Note. — First and second harmonic fit values per declination 
for the second energy band. 




Fig. 6. — Figure (a) shows the one-dimensional projection in 
right ascension a of the first energy band (20 TeY) of two- 
dimensional cosmic ray map in Figure 5-a. Figure (b) shows the 
one-dimensional projection in right ascension a of the second en- 
ergy band (400 TeV) of two-dimensional cosmic ray map in Fig- 
ure 5-b. The data are shown with statistical uncertainties, and the 
black line corresponds to the first and second harmonic fit to the 
data. 
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TABLE 3 

In this table a summary of the sidereal anisotropy energy dependence is displayed. The first column is the median energy 

OF THE COSMIC RAY PRIMARY PARTICLES FOR THE FIRST AND SECOND ENERGY BAND. THE SECOND COLUMN IS THE NUMBER OF EVENTS USED 
IN THE ONE-DIMENSIONAL PROJECTION FROM DECLINATION -24 TO DECLINATION -72. THE VALUES OF THE FIRST AND SECOND HARMONIC 
FITS AMPLITUDES AND PHASES TOGETHER WITH THEIR STATISTICAL AND SYSTEMATIC UNCERTAINTIES ARE DISPLAYED IN COLUMN THREE 
THROUGH SIX. THE LAST COLUMN IS THE \ 2 /ndof FOR THE FIRST AND SECOND HARMONIC FIT TO THE ONE-DII IENSIONAL PROJECTION. 


E Median events A\ SI p <f>\ S1D A 2 SI p <f>2 $ ID X 2 Jrido} 

(TeV) (10 9 ) (IQ""*) (degree) (10~*) (degree) 

20 17.9 7.9±0.1«tot. ± 0.3 S yst. §b".5 ± l-0 s * a t dh 1-lsyst. 5". 9" ± 0.1 a toi. ± 0.4 svs t. 299.5 ± 1.3 s tat i 95/19 

400 0.5 3.7 ± 0.7 stat. ± 0 .Islet. 239.2 ± 10 S 3t at ± 2.7 ± 0.7 stat. ± 0.6 S y 9t . 152.7 ± 7.0 s tat ± 4.2,^. 34/19 


applied to the significance skymaps to improve the sen- 
sitivity to large features. The smoothing search applied 
in this analysis is from 1 to 30 degrees. After smoothing 
is optimized, the significance is then calculated using the 
method of Li & Ma (1983). 

The maximum significant features in the 20 TeV map 
with a 30 degree smoothing are found with an excess at 
(a = 80.8°, S =-49.7°) with a significance value of 40a, 
and a deficit at (a = 219. 7 °, 8 —52.0°) with a significance 
value of -53.5a. Moreover, for the 400 TeV map, two re- 
gions were identified to be significant. The first region is 
an excess at (a = 256.6 ° ,5 —25.9°) with a significance 
of 5.3a and an optimized smoothing of 29 degrees, and 
the second region is a deficit at (a = 73.1°, <5 =-25.3°) 
with a significance of -8.6a and an optimized smoothing 
of 21 degree. These significance values do not account 
for the scan for the peak significance in all pixels of the 
sky or the scan over smoothing radii applied to obtain 
an optimal sensitivity to the observed features. We con- 
servatively estimate a trial factor by assuming that all 
scans give statistically independent results. After cor- 
recting for the trials, only the deficit remained significant 
beyond the 5a level, with a post-trial significance value 
of -6.3a. This is the first significant observation of an 
anisotropy in the Southern sky at 400 TeV. The impli- 
cations of this observation is explored in the conclusion 
and discussion sections of the paper. 



Fig. 7. — Figure fa) shows the significance map for the 20 TeV 
energy band plotted v/ith 30 degree smoothing. Figure (b) shows 
the significance map for the 400 TeV energy band plotted with 20 
degree smoothing. 


3.2. Solar Dipole Anisotropy 


Currently there is no detailed theoretical model that 
predicts the observed sidereal anisotropy in the cos- 
mic ray arrival direction distribution. Except for test- 
ing the stability of the Observatory and its time cov- 
erage (see sec. 4), the only effective way to have 
an absolute calibration of the experimental sensitiv- 
ity for the detection of the sidereal directional asym- 
metries is to measure the solar anisotropy from the 
Earth’s orbital motion around the Sun. The so- 
lar anisotropy is well understood and was first re- 
ported in 1986 by Cutler & Groom (1986) and then 
later observed by experiments in the multi-TeV energy 
range (Tibet AS7 (Amenomori et al. 2004, 2006), Mi- 
lagro (Abdo et al. 2009) and EAS-TOP (Aglietta et al. 
2009)). The observed solar anisotropy consists of a dipole 
that describes an apparent excess of cosmic rays in the 
direction of Earth’s motion around the Sun and a deficit 
in the opposite direction. The relative intensity of the 
solar dipole is expressed as 

Jfj = (7 + 2) ^ cos(0 v ), (3) 

where I is the intensity, 7 the differential cosmic ray spec- 
tral index, v the Earth’s velocity, c the speed of light, and 
Q.j the angle between the reconstructed arrival direction 
of the cosmic rays and the direction of motion of the 
observer (Compton & Getting 1935; Gleeson & Axford 
1968). The actual amplitude of the observed solar dipole 
depends on the geographical latitude of the observer and 
on the angular distribution of the detected cosmic ray 
events at the observatory. 

Due to the location of IceCube at the South Pole, the 
sky is fully visible at any given time. Therefore, the solar 
anisotropy is observed in a reference system where the 
location of the Sun is fixed, where the latitude coordinate 
is the declination and the longitude is defined as right 
ascension difference of the cosmic ray arrival direction 
from the right ascension of the Sun (a - a aun ). In this 
reference frame the dipole excess is expected to be at 
270° and the deficit at 90°. 

Figure 8 and Figure 9 show the maps of the cosmic 
ray arrival direction in solar reference frame, for both 
energy samples (20 and 400 TeV) along v;ith with their 
projection onto right ascension relative to the Sun. The 
color scale is the relative intensity value for each pixel 
normalized to unity for each declination band. A fit to 
the projection of relative intensity distribution vs. (a — 
OLsvn) was done using the first harmonic term of eq. 2. 
Table 4 shows the results of the first harmonic amplitude 
and phase along with x 2 jndof of the fit. 

To verify that the experimental observation of the so- 
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TABLE 4 

First harmonic fit values of the solar dipole anisotropy 

TOGETHER WITH THEIR STATISTICAL UNCERTAINTIES FOR THE 
ENERGY BANDS CENTERED AT 20 TeV AND 400 TeV. 


^Median 

(TeV) 

A Isqi 

(10- 4 ) 

<f>l S OL 

(degree) 

X'/ndof 

20 

1.9 i: O.latat. 

267.1 i 3.6jtat. 

53751 

400 

2.9 ± 0.7 s «at. 

272.1 ± 

12/41 


r-i r cm 

■ ICO T3 

VI roof 

— 1 10005 
3€f ,-H« 

0 9935 
I j 0.909 
0.9985 
0Q9D 

(a) 



(b) 


Fig. 8 — Figure (a) shows the IceCube cosmic ray map of the first 
energy band (median energy of 20 TeV) for the relative intensity 
in right ascension from the sun (a — a 8un ). Figure (b) shows the 
IceCube cosmic ray map of the second energy band (median energy 
of 400 TeV) for the relative intensity in right ascension from the 
SUn (a — Qfstin)- 

lar dipole is consistent with expectation, the predicted 
projection of the solar anisotropy is calculated for the 
IceCube location. The expectation of the solar dipole 
was calculated by computing the relative intensity of the 
solar dipole through the cosmic ray plasma (eq. 3). In- 
stead of counting the number of events within a given 
bin in right ascension from the Sun, for each event, after 
time scrambling the data we calculated a mean weight 
corresponding to the expected relative intensity of the 
solar dipole . 

The uncertainties in the cosmic ray spectral index, in 
the reconstructed arrival direction of the events, and the 
spread in the Earth’s velocity over a year were included 
in the calculation of the uncertainty of the expectation. 
The mean spectral index was evaluated using the all- 
particle cosmic ray spectrum from Horandel (2003) and 
the spec;ral index was found to be (7) = 2.67± 0.19. 
The value used for Earth’s velocity was v = 29. S± 0.5 
km/s Williams 2004), where the error takes into ac- 
count the spread between the maximum and the min- 
imum along the elliptical orbit. The angle 0 V between 
the reconstructed direction of the muon events and the 
Earth’s velocity vector at the time the event was detected 
was evaluated by accounting for the experimental point 
spread function. The expected solar dipole distribution, 
including the 68% spread in the uncertainty of the ex- 
pectation, is shown as a shaded band in Figure 9 for the 
low and high energy samples, respectively. The figures 
show t.iat the observations are consistent with the ex- 
pectation in both amplitude and phase for both low and 




Fig. 9. — Figure (a) shows the one-dimensional projection in 
right ascension from the sun (a — a 3un ) of the first energy band 
20 TeV) of two-dimensional cosmic ray map in Figure 8-a. Figure 
b) shows the one-dimensional projection in right ascension from 
the sun (a — ot 3 un) of the second energy band (400 TeV) of two- 
dimensional cosmic ray map in Figure 8-b. The data are shown 
with statistical uncertainties, and the black line corresponds to the 
first and second harmonic fit to the data. 

high energy distribution. This demonstrates the reliabil- 
ity of the analysis to identify anisotropies at the level of a 
few 10~ 4 , which supports the observations of the sidereal 
anisotropy. 

4. SYSTEMATIC UNCERTAINTIES OF THE SIDEREAL 
ANISOTROPY 

In order to assess and quantify the systematic uncer- 
tainties in the sidereal anisotropy for the low and high 
energy samples of the cosmic ray arrival direction dis- 
tribution, we performed two different studies, similar 
to (Abbasi et al. 2010a). First of all, we estimated the 
statistical stability of the result and verified that the ob- 
servation is unaffected by the particular choice of the 
data sample. Then we estimated the possible distortion 
effect on the sidereal anisotropy distribution derived from 
a possible annual modulation of the amplitude of the so- 
lar anisotropy. 

4.1. Data Stability 

To assess the stability of the sidereal anisotropy, checks 
were applied by dividing the full data sample used in this 
analysis into series of two exclusive data sets by splitting 
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both high and low energy data samples in halves based 
on different criteria. A full analysis was done with each 
dataset and the relative intensity distribution in right 
ascension was determined for each of them, along with 
a fit to the first and second harmonic term of eq. 2, and 
compared to the ones from the complete low and high 
energy samples, respectively. 

To check if the anisotropy had a seasonal dependence 
the data were divided into austral summer and austral 
winter sets. The summer set included events collected 
from December to May while the winter set included 
events collected from June to November. Since each 
dataset used in this test did not cover the full year, the 
sidereal anisotropy distribution was contaminated by the 
un-co:npensated solar dipole (see section 3.2). This spu- 
rious effect was accounted for by determining what the 
solar dipole should look like in a sidereal reference frame 
within the two seasonal time periods. In order to do 
so a numerical calculation was performed where, every 
100/is, an event was generated with a unique UTC time, 
and with right ascension from the Sun sampled from the 
all-year experimental solar dipole distributions for each 
energy sample as shown in Figure 9. The corresponding 
distributions in the sidereal reference frame were then 
calculated and subtracted from the observed sidereal dis- 
tribution in each seasonal time interval and the corrected 
sidereal distributions were then obtained. 

To ensure that the sidereal anisotropy was not affected 
by uniform variations in rate, the daily median rate was 
determined and two data sets were selected. One dataset 
containing sub-runs with event rate above the median 
daily value, where a sub-run corresponded to approxi- 
mately 2 minutes of observations, and one with event 
rate below the median daily rate. Once more the anal- 
ysis was then applied to each dataset and the sidereal 
anisotropy distributions for these data sets were deter- 
mined. 

To check whether the measurement is stable against 
the choice of the particular event sample selection, two 
separate sub-run selection tests were applied. The first 
test was done by dividing the sub-runs randomly for each 
day in two halves, and the second by dividing in even- 
and ocd-numbered sub-runs. The arrival direction dis 1 
tribution in sidereal reference frame was then determined 
for each of these data sets. 

For each day good quality runs were selected that sat- 
isfied fundamental data integrity requirements. This run 
selection, along with sporadic data acquisition downtime 
resulted in data collection time gaps which represented 
about 10% of the livetime in IceCube-59. To verify that 
the non-uniform time coverage due to gaps in the data 
was correctly handled by the azimuthal re- weighting pro- 
cedure. a complete analysis was performed on the sub- 
sample of days with maximal data collection time (i.e. 
~24 hr). There were 214 such days during one calendar 
year of IceCube-59 phvsics run. The relative intensity 
was then determined for the cosmic ray arrival direction 
in sidereal reference frame for this data set. 

The sidereal distributions of relative intensity in the 
cosmic ray arrival direction for the low and high energy 
samples and for each of the .above mentioned tests, were 
used to evaluate the spread in the experimental observa- 
tion from the full-year event samples. The gray bands in 
Figure 10 and 11 describe the maximal spread obtained 


from the result of all the stability checks described in this 
section. 



Fig. 10. — The one-dimensional projection in sidereal time frame 
of the two-dimensional cosmic ray map in Figure 6-a for the 20 
TeV band. The data are shown with statistical uncertainties, and 
the black line corresponds to the first and second harmonic fit to 
the data. The gray band indicates the maximal spread from the 
stability checks. 



Fig. 11.— The one- dimensional projection in sidereal time frame 
of the two-dimensional cosmic ray map in Figure 6-b for the 400 
TeV band. The data are shown with statistical uncertainties, and 
the black line corresponds to a the first and second harmonic fit to 
the data. The gray band indicates the maximal spread from the 
stability checks. 


4.2. anti- sidereal time 

The sidereal anisotropy will be distorted by the solar 
dipole unless data are collected within an integer num- 
ber of full years. While the sidereal reference frame is 
defined where the celestial sky is fixed, the solar refer- 
ence frame is defined where the Sun is fixed. This means 
that a sidereal day is on average 4 minutes shorter than 
a solar day, and therefore, while the solar time reference 
frame includes 365.25 days/year, the sidereal time ref- 
erence frame is composed of 366.25 days/year. A static 
point in the solar reference frame will move across the 
sidereal frame and return to the same position on the 
sky in one full year. As a consequence any static solar 
distribution averages to zero in sidereal reference frame 
after one year. 
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The situation however changes if for some reason the 
measured solar anisotropy has, for instance, an annual 
modulation of its amplitude. Since a non-static signal 
in solar reference frame does not average to zero in side- 
real fiame after one vear, particular care is needed to 
account for this possible source of bias in the sidereal 
anisotropy. This introduces a bias in the reference frame 
where one day is 4 minutes shorter than a solar day (i.e. 
the sidereal frame) and an equivalent bias in the refer- 
ence frame where one day is 4 minutes longer than a 
solar day. This defines the so-called anti-sidereal time, 
i.e. a non-physical reference frame obtained by reversing 
the sign of the transformation from solar time to side- 
real time, where the anti-sidereal year consists of 364.25 
days (Nagashima et ai. 1983). The antisidereal reference 
frame can, therefore, be used to quantify the distor- 
tion induced in the sidereal anisotropy (Farley &; Storey 
1954). 

Figure 12 shows the relative intensity of cosmic rays 
arrival distribution in anti-sidereal reference frame (for 
the low and high energy samples). The anti-sidereal 
anisotropy is measured by using a coordinate system 
where the longitude coordinate is defined using the anti- 
sidereal time (<*a,s)- The figure also shows a fit to the 
observed distributions with the dipole term of eq. 2 and 
Table 5 shows the fit results. Both the low and high 
energy samples show no significant observed amplitude 
in the anti-sidereal time. The uncertainty in the first 
harmonic amplitude and phase derived by the study of 
the anti-sidereal distribution was found to be within the 
statistical and systematic errors determined from the sta- 
bility tests. 

The results of all the systematic checks described in 
section 4.1, along with the estimate of the distortion in 
the sidereal anisotropy, based on the anti-sidereal dis- 
tribution, were collectively used to estimate the global 
systematic uncertainties in the sidereal anisotropy fit pa- 
rameters. Adding these systematic uncertainties the first 
and second harmonic amplitude and phase of the side- 
real anisotropy for the low and high energy samples are 
summarized in Table 3. 


TABLE 5 

First harmonic fit values of the anti-sidereal anisotropy 
FOR THE ENERGY BANDS CENTERED AT 20 TeV AND 400 TeV. 


B Mediin 


&1-ASID 

X 2 /ndof 

(TeV: 

(io-4) 

(degree) 


20 

0.4 ± 0.1 

1.5 ±18.5 

29/21 

400 

0.5 ±0.7 

324.6 ± 75.4 

17/21 


4.3. IceCube -40 String Sidereal Anisotropy 

In addition to the previously discussed systematic 
checks, an important cross-check is applied by looking 
at the result obtained from the previous year using the 
data collected from IceCube in its 40-strings configura- 
tion (IceCube-40) from May-2008 until May-2009. The 
same analysis described in this paper v/as applied to the 
IceCube-40 experimental data, along with the energy 
sample selection described in section 2.2. 

The sidereal anisotropy observed at 20 TeV with 
IceCube-40 is found to be consistent with the reported 




Fig. 12. — Figure (a) shows the projection in olas of the relative 
intensity of cosmic ray arrival distribution using the anti-sidereal 
time for the low energy sample (median energy of the primary 
cosmic ray particle of 20 TeV). Figure (b) shows the projection in 
a as of the relative intensity of cosmic rays arrival distribution for 
the high energy sample (median energy of the primary cosmic ray 
particle of 400 TeV). An anisotropy in the anti-sidereal reference 
frame is related to a distortion of the sidereal anisotropy induced 
by an annual modulation of the solar dipole amplitude. 

observation with IceCube-22 (Abbasi et al. 2010a) and 
with that observed using the IceCube-59 string configu- 
ration. Moreover, the relative intensity distribution for 
IceCube-40 as a function of right ascension for the 400 
TeV band is also consistent v/ith the distribution ob- 
tained with IceCube-59. Figure 13 shows the projection 
in right ascension of the relative intensity distribution 
at primary median energy of 400 TeV for both IceCube- 
40 (in red) and IceCube-59 (in black). The line corre- 
sponds to the first and second harmonic fit to the data 
of IceCube-59 in black and IceCube-40 in red. The gray 
band indicates the estimated maximal systematic uncer- 
tainties of IceCube-59. The results obtained with the two 
detector configurations are consistent within the statis- 
tical and systematic fluctuations. The stability of the 
result over different detector configuration supports the 
conclusion that the anisotropies observed at 20 and 400 
TeV with IceCube-59 are real. 

5. CONCLUSIONS 

In this paper we presented the results on the large scale 
cosmic ray sidereal anisotropy, based on a total of 33 x 10 9 
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Fig. 13.— This figure shows the IceCube-59 and IceCube-40 one- 
dimensional projections in sidereal time in black and red marliers 
respectively at 400 TeV. The data are shown with statistical uncer- 
tainties for error bars. The line corresponds to the first and second 
harmonic fit to the data of IceCube-59 in black and IceCube-40 in 
red. The gray band indicates the estimated maximal systematic 
uncertainties of IceCube-59. 

muon events collected by IceCube-59 from May 2009 to 
May 2010. In particular we showed the relative intensity 
in the arrival direction distribution at primary particle 
median energy of about 20 TeV and 400 TeV as shown 
in Figure 5. 

The relative intensity distributions as a function of 
right ascension is fitted with a sum of first and second 
harmonic terms (eq. 2). The amplitude and phase at 
20 TeV and 400 TeV are summarized in section 4. The 
observation of the sidereal anisotropy in the cosmic ray 
arrival direction is supported by the determination of the 
solar dipole expected from the Earth’s revolution around 
the sun. The observed solar anisotropy agrees in am- 
plitude and phase with the expectation in both energy 
bands. Moreover, the sidereal anisotropy is also sup- 
ported by a number of data stability checks. One of 
these checks consisted of analyzing the data samples in 
the ami-sidereal.time frame where no significant signal is 
observed. The observation of the solar dipole along with 
the absence of a signal in the anti-sidereal time frame in 
addition to all the stability tests, ensure the reliability 
of the sidereal anisotropy measurement for both 20 TeV 
and 400 TeV primary energy event samples. 

The sidereal anisotropy observed at 20 TeV with 
IceCube-59 is consistent with the previously reported 
observation with IceCube-22 (Abbasi et al. 2010a), thus 
providing a confirmation of a continuation of the arrival 
distribution pattern observed in the Northern equato- 
rial hemisphere in the multi-TeV energy range. On the 
other hand the sidereal anisotropy observed at 400 TeV 
shows a significant relative deficit region in right ascen- 
sion, where the excess is observed at median pri- 

mary energy of 20 TeV. In addition, the relative deficit 
region at low energy seems to have disappeared at me- 
dian primary energy of 400 TeV as shown in Figure 7. 
The observed anisotropy at 400 TeV shows substantial 
differences v/ith respect to that observed at 20 TeV. 
Moreover, it does not show a continuation of the obser- 
vations reported at high energies in the Northen hemi- 
sphere (Amenomori et al. 2006), (Aglietta et al. 2009). 
This is the first significant anisotropy observed in cos- 
mic ray arrival distribution in the 400 TeV range in the 


Southern hemisphere. 

The sidereal anisotropy at 400 TeV also appears to be 
present in the data collected during the 40-string Ice- 
Cube physics runs. The persistence of the anisotropy 
in IceCube-40 and IceCube-59 is an important verifica- 
tion that the anisotropies observed are not dependent 
on the detector configuration nor on the period the data 
were collected. Using events collected with the complete 
IceCube observatory (86-strings) will enable us to signifi- 
cantly improve the statistical power in the determination 
of sidereal anisotropy at a few hundreds TeV primary en- 
ergy. 

6. DISCUSSION 

The origin of the sidereal anisotropv is still unknown. 
It is believed that a possible contribution to this ob- 
served anisotropy might be from the Compton- Getting 
effect. The Compton- Getting dipole anisotropy we ex- 
pect to see in this analysis is determined from Monte 
Carlo simulation and should appear with a maximum 
in the one-dimensional projection in right ascension be- 
tween 290°and 340° and a deficit between 110° and 160° 
with an amplitude of ~ 0.13%. In this model the cosmic 
rays are assumed to be at rest with respect to the Galac- 
tic center. The sidereal anisotropy from both energy 
samples do not appear to be consistent with expected 
Compton- Getting model (Compton & Getting 1935) nei- 
ther in amplitude nor in phase. However, it is possible 
that the Galactic cosmic ray rest frame has a smaller 
relative velocity and a different direction with respect to 
the one hypothesized in (Compton &: Getting 1935). The 
cosmic ray rotation with respect to the Galactic center 
is complex and unknown, therefore, in this case we can 
only conclude that the cosmic rays are not at rest with 
respect to the Galactic center. 

It is also worth noting that when describing the Galac- 
tic cosmic ray propagation through diffusion models 
the large scale anisotropy is an important observable. 
The determination of cosmic ray anisotropy at median 
energy of 400 TeV could enable us to obtain an im- 
proved theoretical description of the diffusion processes 
of Galactic cosmic ray’s energy ranges closer to the 
knee (Berezinskii et al. 1990). 

We are continuously analyzing events from IceCube 
with updated configurations. IceCube construction is 
now completed with 86 strings deployed with a volume 
of km 3 in January of 2011. With the higher statisti- 
cal power expected from the observed cosmic ray muons 
we will be able to improve our understanding of the 
anisotropy and its energy dependence closer to the knee 
region. This will further our understanding of the propa- 
gation of cosmic rays and help to eventually reveal their 
sources. 
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